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^ ■ Abstract 

' , We analyze the cascade decays of the scalar quarks and gluinos of the Minimal 

Supersymmetric extension of the Standard Model, which are abundantly produced 
at the Large Hadron Collider, into heavier charginos and neutralinos which then 
decay into the lighter ones and charged Higgs particles, and show that they can 
have substantial branching fractions. The production rates of these Higgs bosons 
can be much larger than those from the direct production mechanisms, in particular 
for intermediate values of the parameter tan/3, and could therefore allow for the 
detection of these particles. We also discuss charged Higgs boson production from 
direct two-body top and bottom squark decays as well as from two- and three-body 
gluino decays. 



1. Introduction 

The most distinctive signature of an extended Higgs sector, compared to the Standard 
Model (SM) where only one scalar doublet is needed to break the electroweak symme- 
try leading to a single neutral Higgs particle, is the discovery of charged Higgs bosons. 
For instance, in the Minimal Supersymmetric (SUSY) extension of the Standard Model 
(MSSM) |p, two Higgs doublets are present, leading to the existence of a quintet of scalar 
particles: two CP-even neutral Higgs bosons h and H, a CP-odd neutral Higgs boson 
A and two charged Higgs particles H 0. SUSY constraints on the Higgs spectrum 
impose that the charged Higgs boson mass is related to the pseudoscalar Higgs mass, 
Mjj± = M\ + My V . With the present experimental limit on the A boson mass, Ma ^ 93.5 
GeV 0, this leads to a bound M H ± > 120 GeV. In fact, in the popular minimal Super- 
gravity models (mSUGRA) with universal boundary conditions at the Grand Unification 
scale and where the electroweak symmetry breaking is induced radiatively j|], the H 
boson, as well as the H and A neutral Higgs particles, tend to be rather heavy, with 
masses of the order of a few hundred GeV 0. The states are therefore kinematically 
accessible only at the LHC 0, at future e + e~ linear colliders [[jj or muon colliders 0. 

The discovery of H bosons at the LHC through the standard processes is rather 
difficult, if not impossible in some areas of the MSSM parameter space ||. This is mainly 
due to the fact that the production rates are controlled by the charged Higgs boson Yukawa 
couplings to up- and down-type fermions. Using the notation of the first generation, the 
latter are given by 0: 

H + [cot (3 m u UidjL + tan (3 m d Uidj R ] , (1) 



V2M, 



w 



where V^ is the CKM matrix with Ui = u, c, t, di = d,s,b and tan/3 = v 2 /vi is the ratio of 
the vacuum expectation values of the two Higgs doublets needed to break the electroweak 
symmetry in the MSSM. For values tan/3 > 1, as is the case in the MSSM, the couplings 
to down-type (up-type) fermions are enhanced (suppressed). Only the couplings to the 
top and bottom isodoublet quarks are therefore important, in particular for small ~ 1 
and large ~ m t /m b values^. 

A light charged Higgs particle can be searched for at the Tevatron in top decays 
through the process pp — > ti, with at least one of the top quarks decaying via t — > H + b, 
leading to a surplus of r leptons due to the H ± —> r ± z/ decay, an apparent breakdown of r- 
\i— e universality. For small and large values of tan /3, the branching ratio BR(t — > H + b) is 



large and would allow for the detection of the signal at the Tevatron |]T0| . The situation for 
intermediate values of tan/3, where the H~tb coupling is small, leading to a tiny t -^ H + b 
branching ratio, is rather difficult and should await for the LHC. Indeed, detailed analyses 
of the ATLAS and CMS collaborations have shown that the entire range of tan /3 values 
should be covered for M H ± < m t using this process H. 



interestingly, these two regions of tan/3 are favored by Yukawa coupling unification. However, the 
experimental bound on the lightest h boson at LEP2, Mh > 113.5 GeV || in the decoupling regime 
where the H^ bosons are heavy, rules out the low tan/3 scenario. 



For charged Higgs bosons with masses M#± > m t , the two production mechanisms 
which potentially have sizeable cross sections at the LHC are [|TT| Q: 

pp - gb(gb) - tH- (tH+) 

PP — ► 99 /qq — ► tH~b + iH + b (2) 

The signal cross section from the 2 — > 2 mechanism (76 — ► £/7~, where the b quark is 
obtained from the proton, is 2-3 times larger than the 2^3 process gg/qq — * tbH~, 
where the H~ boson is radiated from a heavy quark line. When the decays H + — > tb and 
t — > Wb take place, the first process gives rise to 3 6-quarks in the final state while the 
second one gives 4 5-quarks; both processes contribute to the inclusive production where 
at most 3 final 6-quarks are tagged^. However, the cross sections are rather small: even 
for the extreme values tan/3 = 2 and 40, they hardly reach the level of a picobarn for 
a charged Higgs boson mass Mh± = 200 GeV. For intermediate values of tan j3 and/or 
larger H masses, the cross sections are too small for these processes to be useful. For 
instance, for the value tan/? = 10, the cross section is below the level of a few femtobarn 
for M H ± > 250 GeV. 

Other mechanisms for H ± production at hadron colliders are the Drell-Yan type pro- 
cess for pair production through 7 and Z boson exchange, qq — > H + H~ [I3j , the gluon- 



gluon fusion process for pair production, gg — > H + H J14l , and the associated production 
process with W bosons in gg fusion and qq annihilation, qq,gg — > H^W 11 ||15| . However, 
the production rates are rather small at the LHC: for the quark-antiquark processes 
because of the low quark luminosities at high energies and for the gluon-gluon fusion 
processes because they are induced by loops of heavy quarks (and squarks) and are thus 
suppressed by additional coupling factors. The cross sections are at the femtobarn level 
for large enough charged Higgs boson masses, and are therefore too low to be easily useful 
in the complicated and hostile environment of hadron colliders. 

In this paper, we show that there is a potentially large source of the H ± bosons of 
the MSSM at the LHC: the cascade decays of squarks and gluinos, which are abundantly 
produced in pp collisions, thanks to their strong interactions^]. Squarks and gluinos can 
decay into the heavy charginos and neutralinos, xt ■> X3 an d X% an d if enough phase space 
is available, the latter could decay into the lighter charginos/neutralinos, X11X1 an d X% 
and charged Higgs bosonsf]: 

pp^gg,qq,qg -» xhx^Xi + x 

- xf,xlx° 1 + H ± +X (3) 



2 However, the two processes have to be properly combined to avoid double counting of the contribution 
where a gluon gives rise to a bb pair that is collinear to the initial proton. The cross section of the inclusive 
process in this case is mid-way between those of the two production mechanisms eqs. (2) B. 

3 These decays have been discussed in the past |1(| for charged Higgs bosons with masses below ~ 150 
GeV and which can thus also be produced in top quark decays. 

4 We will not specifically consider here the similar production of neutral Higgs particles, $ = h, H, A, 
since the later can be produced abundantly in standard processes such as as gluon-gluon fusion, gg — > $ 
and associated production with heavy quarks gg, qq — * 6M>, tt$ [||. 



These processes are similar to the ones with cascade decays of strongly interacting SUSY 
particles into the next-to-lightest neutralino x® which then decays into the lightest h 
boson and the lightest neutralino [which is expected to be the lightest SUSY particle in 
the MSSM], a process which has been discussed in the literature; see for instance Ref. ]17[ . 
Charged Higgs bosons could also be searched for, if kinematically possible, in the 
direct decays of heavy third generation squarks into their lighter partners, 

Q-^Q'H* with Q,Q' = i,b (4) 

or in direct gluino three-body decays into heavy quarks, their partners squarks and H ± 
bosons, 

~g^Q'QR ± with Q = i,b (5) 

The remainder of the paper is organised as follows. In the next section, we discuss the 
production cross sections of squarks and gluinos at the LHC as well as their main decay 
modes. In section 3, we analyse the decays of the heavy charginos and neutralinos into 
the lighter charginos/neutralinos and Higgs or gauge bosons and estimate the production 
cross sections for H ± bosons at the LHC in some specific scenarii. In section 4, we discuss 
the direct production of H ± bosons in the two-body decays of top and bottom squarks 
and three-body decays of gluinos. A short conclusion is given in section 5. 

2. Production and decay modes of Squarks and Gluinos 

2.1 Production cross sections 

In proton-proton collisions, gluino pairs are produced through qq annihilation and gluon- 
gluon fusion, qq, gg — > gg. Squark pairs can be produced through t-channel exchange, 
qq — > qq, while squark-antisquark pairs are produced in both qq annihilation and gluon- 
gluon fusion. Finally, mixed squark-gluino production proceeds through s- and t- channel 
gq annihilation. The production cross sections at the tree-level are given in Ref. JR 



They have however, to be supplemented by the next-to-leading order QCD radiative 
corrections, which stabilize the theoretical predictions. The i^-factors can be rather 



large, ranging from K = onlo/olo ~ 1 to 2 [[19], depending on the ratio of the squark to 



gluino masses. Taking for the renormalisation and factorization scale the average mass m 
of the two produced sparticles results in a conservative estimate of the total production 
cross section. We will therefore not include the i^-factors in our analysis and set the scale 
at which the cross section is evaluated to be the average mass of the final particles. 

The cross sections for pair production and associated production of squarks and gluinos 
are shown in Fig. 1 as functions of the squark and gluino masses. The CTEQ3L p0 | 



parameterization of the parton densities has been used. In Fig. la, we choose m g = rriq 
while in Figs, lb and lc, we take respectively, rriq = 1.2m g and m g = 1.2rriq. We assume 
the left- and right-handed scalar partners of the light quarks, including the 6-quark, to 
be degenerate in mass and we sum the individual cross sections. We show separately 
in Fig. la, the total cross section for the pair production of the lightest top squark, 
a qqi99 ~~ * tit*)) which depends only on m^ at the tree-level. 
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Figure 1: Pair and associated production cross sections at the lowest order for the strongly 
interacting SUSY particles at LHC with y^s = 14 TeV as a function of the final state 
masses for representative choices of squark and gluino masses, a^ and <jq g include the 
contributions from a^* and Uqg respectively; a^j* is presented only in Fig. la. 



We see that the largest cross section is due to squark and gluino associated production, 
in which squarks and anti-squarks of the five light flavors have been added. It ranges 
from a— ~ 50 pb for a gluino mass nig ~ 500 GeV to O (1 pb) for m g ~ 1 TeV. In the 
low mass range, it is followed by the cross section for gluino (squark) pair production if 
the gluino is lighter (heavier) than squarks. The total cross sections decrease quickly with 
increasing final state masses. In the case of top squarks, the total production cross section 
is an order of magnitude smaller than the one for squark-antisquark pairs, a^j* ~ j2 a qq* 
if all squarks, including t\, have the same mass; the cross section for the production of t 2 , 
which is expected to be heavier than t\ because of the generally large mixing in the stop 
sector, is smaller than cr^j*. 

Summing up all cross sections, one obtains a q+ g ~ 110 (3) pb for m~ g ~ nig ~ 500 
(1000) GeV. This means that with the expected integrated luminosity at the LHC, J r £ ~ 
300 fb _1 , a total of 3 • 10 7 to 10 5 events can be collected by the ATLAS and CMS 
collaborations in a course of few years. It is therefore tempting to use this very large 
sample of events to look for particularly interesting decay channels of squarks and gluinos 
which, even if they have branching ratios below the percent level, would still lead to a 
rather large number of final state events that could be studied in detail. 

2.2 Squark Decays 

The main decay modes of squarks, if they are lighter than the gluino, will be into their 
partner quarks and neutralinos, q~i — > qx°j [7=1-4], as well as quarks and charginos, 
Qi ~^ Q'xt b = l~2]. Taking into account the mass of the final quark, which would be 
appropriate in the case of top squark decays, the partial decay widths are given at the 
tree-level by [the QCD corrections to these decay modes have been calculated in Ref. pi 



aA5(/i2,/i2 ) 
r(ft -► QXj) = 4 —mqi 

aA5(^„^2 + ) 

r(ft -► q'Xj ) = : —m. 



(rfj + % )(1 — A*? — /4°) ~ ^Wii^xfxi 



>i, 



a \j + b ij )(1 - /4 ~ t X l+) - Aa 1j b ij^q'^ 



(6) 



where X(x, y) = l + x 2 + y 2 — 2(xy + x + y) is the usual two-body phase space function with 
the reduced masses \ix = fn x /mg. and e Xj is the sign of the eigenvalue of the neutralino 
X°. In terms of s 2 ^ = 1 — c^ = sin 2 Wl the squark electric charge, weak isospin and the 
mixing angle 9 q which turns the left- and right-handed states into the mass eigenstates 
<?i = c e q QL + se q qR and 52 = —se q qL + C9 q qn, one has for the couplings among neutralinos, 
quarks and squarks: 



V2M w s w \ c e , 
b q j2 J V2M w s w \ so, 







(7) 



with r u = Zji/ sin /3 and r d = Zj 3 / cos (3 for up and down-type fermions, and 
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(8) 



while for the couplings among charginos, fermions and sfermions, qi — q' — Xji one nas f° r 
up-type and down-type sfermions: 
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(9) 



(10) 



In these expressions, U, V and Z are the (real) diagonalizing matrices for the chargino 
and neutralino states J22] with: 



Z ix — ZuCy/ + Z i2 S\Y , Z i2 — —ZiiSy/ + Z i2 cw , Z a — Zj 3 , Z iA — Z,i 



/•i 



(11) 



If squarks are heavier than the gluino, they can also decay into gluino-quark final states, 
for which the partial decay width is given by: 



r(& -> ?<?) 



2a s A3(//2 j// |) 



m,, 



1 — A*2 — A*l — ^Uf-JXgfXg 



with the same notation as previously and the squark-quark-gluino coupling are: 



''!</ 



6| g = sin 6 q , og g 



-6f^ = cos e q 



(12) 



(13) 



We will now discuss some scenarii for these decay modes, starting with the case of the 
scalar partners of light quarks and continuing with the special case of top squarks. 

2.2.1 The case of the scalar partners of light quarks 

If squarks are heavier than the gluino, they will decay most of the time into quark plus 
gluino final states. This is essentially due to the fact that these are strong interaction 
decays, compared to their weak interaction decays into charginos and neutralinos. The 
large value of the strong coupling constant, a s /a ~ 10, makes the decays into gluinos 
an order of magnitude larger. In the opposite case, rriq < m g , the right-handed squarks 
will decay [for small quark masses] only into quarks and neutralinos, while left-handed 
squarks decay into both charginos and neutralinos. Two scenarii are possible, depending 
on the chargino/neutralino textures: 



7 



(i) Gaugino-limit: If the lighter chargino and neutralinos are gaugino-like, that is if 
the higgsino mass parameter is much larger than the wino and bino mass parameters, 
|/i| ^> M 1)2 with Mi and M 2 related by the GUT constraint M 2 ~ 2 Mi, the masses are 
such that m x o ~ m ± ~ 2m x o ~ M 2 while the heavy chargino and neutralinos have 
masses, m x o ~ m x o ~ m ± ~ |/i|. Squarks will have then the tendency to decay into the 
lighter ino states not only because of a more favorable phase space, but also because for 
the partners of the light quarks, the higgsino component of the quark-squark-ino coupling 
[which is proportional to m q , see the previous formulae for the couplings] is very small. 
Therefore squarks will decay dominantly into the lighter charginos and neutralinos. 

(ii) Higgsino-limit: If, the lighter chargino and neutralinos are higgsino-like, that is 
Mi o ^> I Ail, the trend is reversed and one would have the mass hierarchies, m v o ~ m + ~ 
m x o ~ |/i| and 2m x o ~ m ± ~ m x o ~ M 2 . If allowed by phase space, i.e. for trig > M 2 , 
(left-handed) squarks will decay into all possible neutralino and chargino combinations 
in principle. However, because the higgsinos couple proportionally to the quark masses, 
the partial decay widths into the lighter chargino and neutralinos are tiny, and squarks 
will dominantly decay into the heavier chargino and neutralinos, q — > q'x 2 ? QX3 an d qxl- 
[In the mixed region, M 2 ~ |/z| decays into all charginos and neutralinos are possible if 
phase space allowed.] 

The two scenarii (i) and (ii) are exemplified in the upper and lower panels of Fig. 2, re- 
spectively, where we display the "average" squark branching ratios into gluinos, charginos 
and neutralinos as a function of the squark mass for the value tan/5 = 10. This "average" 
branching ratio is defined as the branching fraction of the decay of any squark [except 
for top squarks] into a given neutralino or chargino [or gluino] final state; this means 
that we sum all possibilities for left- and right-handed squarks as well as for up-type 
and down-type squarks keeping proper track of flavors and chirality to reach a given final 
state. [Note that here we neglect the small effect of the 6-quark Yukawa coupling which 
is not enough enhanced for the value of tan/3 we are using, and treat the b squark in 
the same footing as the first and second generation squarks]. The gluino mass is fixed to 
nig = 3M 2 , which is typically the case in models with unified gaugino masses at the GUT 
scale which lead to the tree-level relation m g ~ M 3 ~ 3M 2 at the low-energy scale. 

We see that in scenario (i), where the wino and higgsino mass parameters are set to 
M 2 = 150 GeV and /i = 400 GeV, squarks will decay mainly into the lighter chargino Xi 
and neutralinos Xi 2 m ^ ne l° w mass range, i.e. trig < 450 GeV, with branching ratios of 
the order of 50% for the charged and neutral decay channels, while the strong interaction 
q — > qg decay channel becomes by far dominating above the gluino mass threshold, with a 
branching ratio reaching 90% at large squark masses. In scenario (ii), where the higgsino 
and wino mass parameters are set to yU = 150 GeV and M 2 = 400 GeV, the gluino is much 
heavier than the squarks, m g > 1.2 TeV, and squarks will decay almost 100% of the time 
into the heavier chargino and neutralinos, with a dominance at large squark masses, of 
the charged decay mode, q — > q'x 2 > which reaches a branching fraction of 50% due to the 
usual dominance of the charged currents over the neutral currents. 

Thus, there exist situations in which heavy squarks can decay into the heavier chargino 
and neutralinos with significant rates. 
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Figure 2: Branching ratios of squarks decaying into gluinos, charginos and neutralinos as 
a function of their mass for tan/5 = 10 and rrig = 450 GeV (top) 1.2 TeV (bottom). The 
gaugino mass is fixed to M2 = mg/3 while the value of the parameter \x is 400 GeV and 
150 GeV in scenarii (top) and (bottom) respectively. 



2.2.2 The case of the top squarks 

The case of top squarks is special: because of the large value of top quark Yukawa coupling, 
there is a sizeable splitting between the two states t 2 and ti, and the latter is in general 
much lighter than all other squarks. In some cases, gluinos can be lighter than the scalar 
partners of light squarks, but heavier than t\. In particular, one can have the mass 
hierarchy rrig > m g > mi + m t and the squarks q will decay almost exclusively into 
gluinos and quarks as discussed previously, and the gluinos will decay into the only two- 
body decay mode which is allowed, i.e. g — > tt\. This means that all strongly interacting 
particles produced at the LHC could decay dominantly into top and stop final states. 

The branching ratios for t\ decays into x + b an d X°t final states are shown in Fig. 3 
for a mass of raj = 600 GeV and tan/? = 10. In the upper (lower) panel, they are shown 
as functions of /z(M 2 ) for a fixed value of M 2 = 150 GeV (/i = 400 GeV). As can be 
seen, there are large regions where the decay t\ —* b\t is dominating. In particular, one 
can see that for small values of /i (~ 150 GeV) and large values of M 2 (~ 400 GeV), 
i.e. in the higgsino region, BR(ti — > b\ 2 ) can reach the level of 50%. But in contrast to 
the scalar partners of light quarks, the branching ratio can also be large in the gaugino 
region, M 2 = 150 GeV and fi = 400 GeV, since now the couplings of top/stop states to 
higgsinos are enhanced by m t . The decays into neutralinos, t\ — > txg 4 are also sizeable 
in the gaugino region, but the phase space is limited because of the large value of m t . 

2.3 Gluino decays 

If gluinos are heavier than squarks, their only relevant decay channel will be into quark 
plus squark final states, g — > qq. The partial decay width, including the quark mass and 
the squark mixing angle to take into account the possibility of decays into top squarks, 
using the notation introduced previously, is given by [here, nx = mx/mg]' 



rO^gft) = m g 



1 +Va- V& + ^aUf-Ugflg 



(14) 



In the case, where the gluino is lighter than the squarks, it will mainly decay into final 
states involving a quark-antiquark pair and charginos or neutralinos. The Dalitz density 
for the decay, taking into account the masses of the final fermions, is given by [E3| : 



(LI ^ 7 """ ^ drfc+dit+dn, 



dx u dxd w k 327T 



-(g^Xjud) = 1 ^ z m~ g ^ 



kl=\ 



■^d i JT^ud 



(15) 



dri = 



kl 



x (1 + fi x + fi d + fi u ) - (fi x + fi d )(l + fi u )) - 2y/jr a {a d kg b d r g + af~ g b d k§ )(fi d + [i x - «)] 
+ 2{a%b% + a%b%)^jr x -[T d [{a d k ~ga d g + b%b%){u - 1 - ft u ) - 2 v ^(a d kg b d g + a d g b d kg )} 
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Figure 3: Branching ratios for top squark decaying into charginos and neutralinos as a 
function of fj, with M 2 = 150 GeV (top) and M 2 with fj, = 400 GeV (bottom). We have 
fixed the other parameters to tan/5 = 10 and m^ = 600 GeV. 
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dr 1* = ( - - ,f, ( - - 7 ft (« + « [(«««« + « (-*■ + * 

I— A*d — Z^iifc + t){—Hd — Hui + cj I 



x (1 + // x + /i u + fi d ) - (fi x + /i„)(l + fi d )) - 2y/jfi(a% s bY s + a^ g blg)(n u + /x x - £)] 






.« „d z,d i j,u lm iA „d ' 



(tt + 1 - fiu - fM d ) + Kga] k af s b% + b u k - g b] k bf- g a^}^(t -n x - /i M ) 

+ [a"*4 a ?i 6 « + %^°WV^5(fi - /x„ - 1) - 2[a>^fcJ + afca^ftj] 

v^v 7 ^^ + [tfkbtgbfgtfi + a u kg a? g a d 3l b] k }^(u - /i x - fx d ) 

+ K4 9 KA + <~A b lkbt 9 ](ui -» x - w d ) + [al g a) k a%a% + b u kg b u ]k b d g b%}^jr x 

(u + i - ii x - 1) + [aj^ajij + ^ fe &%]VW# - ^ - 1)} (16) 

where x u = 2E u /m g , x<* = 2E d /m g are the reduced energies of the final quarks, \xx = 



m 



/m~ the reduced masses and t = (p g — p u ) 2 /ml = 1 — x u + /i u , it = (p g — pd) 2 ! /m| 



1 — Xd + /id. The squark-quark couplings to charginos and neutralinos, a^ and 6^ and the 



couplings to gluinos a q ~ and b q ~ have been given previously. The fully integrated partial 
decay width, in the case of massless final state quarks, can be found in Ref. |24| . 



The various gluino branching fractions are shown in Fig. 4 as a function of /x for a 
gluino mass m g = 800 GeV and for squark masses fixed to rrig = 1.2771^ = 1.2m g . Two 
different scenarii are exemplified: in the upper panel of Fig. 4, the approximate GUT 
relation M 2 = m g /3 is adopted while in the lower part of Fig. 4, this GUT relation is 
relaxed and we set M 2 = 150 GeV, i.e. with gaugino-like Xi>X?2 f° r l ar g e A 4 values. 
[We however, keep the equality of the wino and bino mass parameters at the GUT scale, 
leading to M 2 — 2M\ at low energies]. In both scenarii, the top squark is lighter compared 
to other squarks [which is generally the case, as discussed previously], m^ = m g , and its 
virtuality in the decay of the gluinos is smaller. Since, in addition, the top squark has 
large couplings to the higgsino states, the decay modes g —> tbxi and 9 -^ ttxl [and 
the charge conjugate states] are more important, the former being in general dominant 
because of the larger phase-space and the stronger charged current couplings. 

In the universal scenario with M 2 ~ m g/3 ~ 270 GeV, the branching ratio for the 
decay g — > tbx 2 is at the level of ~ 30% for small values of \i < 300 GeV, where x\,i 
are mixtures of gauginos and higgsinos, while the branching ratios of the decays into top 
quarks and X34 final states are smaller by a factor of 3. The branching ratios decrease 
with increasing /j, since xt 4 and X2 become heavier and the phase space for the decay is 
reduced. However, even for values of \x around 400 GeV, BR(g — > tbx 2 ) stays above the 
level of 10%. In the non-universal scenario with M 2 = 150 GeV, xf are almost higgsino- 
like for relatively large /j values, /j > 300 GeV. BR(g — > btx 2 ) is larger than what was 
previously for small /j values while for large \x values, it is similar to the previous case. 
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Figure 4: Branching ratios of gluinos into electroweak gauginos as functions of ii for 
tan/3 = 10. The gluino mass is fixed to m g = 800 GeV while the squark masses are 
rriq = 1.2'mg = \.2m~ g . The wino mass parameter is fixed to M2 = 2M\ = m g /3 (upper 
curve) and M2 = 2M\ = 150 GeV (lower curve). 
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3. Decays into Charged Higgs Bosons 



3.1 Two— body decays of heavy charginos and neutralinos 

The heavier chargino and neutralinos will mainly decay into the lighter chargino and 
neutralino states and, if enough phase space is available, into sfermion-fermion pairs. 
The partial decay widths of the two-body decays, [including the possibility of massive 
fermions in the last case], are given by: 



r(xi - ffi 



aN r . 



m, 



((4) 2 + ( & i) 2 ) C 1 - »h + av) + 4 \^f a ij b ij A ^(/ i /' h) ( 17 ) 



r(Xz -» XjV) = -m Xi \3(n x .,n v ) {-12^/JI^GfiyG 



R 

jiV 



+ 



{G% v f + (Gf iV ) 2 ] (1 + ^ - ii V ) + (1 - A** + /iy)(l - \x Xi - \i v )\if\ (18) 



T( X i -> Xi^fe) = g m Xi A2 (n Xj ,HH k ) { [(Gf jk ) 2 + {Gf jk f\ (1 + /^. - A^J 

+ 4 v /^Gy fc Gy fc | (19) 

The couplings among charginos, neutralinos and fermions/sfermions have been given pre- 
viously, while the chargino and neutralino couplings to the Higgs and the gauge bosons 
are given by [H k — h,H,A and H ± for k — 1, 2, 3, 4]: 



G 



G% 



L,R 



L,R 



with G * jW ~ vt^[~ ZuV i' 2 + ^ Z ^ V ji] 



(20) 



n L,R 

x t xjz 



r L,R 



with 



G L - = 1 
G 11 = 1 



-\Vi2~Vj2 ~ VnVji + SijSfy, 



2~Ui2Uj2 — UnUji + OijS 



w 



(21) 



G x°x°z- G ^ wltn £* 



1 



■[^3^. 



ZnZj4\ 



2swcw L^ w ^i3 

+ o„ [ZizZjz — Z iA Zji\ 



2sw c w ' 



(22) 



G 



L,R 

X°iXt H+ 



G%* with f ^ 



T ij4 



ZjaVh + -^ (Zj 2 + t&n9 w Zji) Vi2 
Zj?,Uii — -^ (Zj2 + t&nOwZji) Ui2 



G 



L,R 

~*-i Xj k 



G ijk Wlth 



Gi i k ~~ 731^ [ e fc^?'i^2 - d k Vj2Un] 
G ^ = 7^ [ekVnU j2 - 4 W,i] e fc 



(23) 



G x°x?^ fe ~ ^w* Wlth r7 R , = -^ 



^ufc = l^F (^' 2 - tan 9 w Zji) (e k Z i3 + 4^) + * <-> j 
u'fc ~~ 2^7 (^' 2 ~~ tan 9wZji) (e k Z i3 + <4Z i4 ) e fc + i <-> j 



(24) 
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where ei i2 = — e 3 = 1 and the coefficients e^ and dk read 

ejdi = c a / - s a , e 2 /d 2 = -s a / - c a , e 3 /d 3 = sp/cp (25) 

The branching ratios of the heavier chargino xf an d X34 m ^° the lighter ones xt 
and Xi 2 an d gauge and Higgs bosons are shown in Figs. 5 and 6 for, respectively, the 
two scenarii (i) gaugino-limit and (ii) higgsino-limit discussed previously. Squarks and 
sleptons are assumed to be too heavy to play a role hereQ. To analyze them, it is useful 
to discuss first the couplings of charginos and neutralinos to the Higgs and gauge bosons. 

The Higgs bosons couples preferentially to mixtures of gauginos and higgsinos. This 
means that in the gaugino-like or higgsino like regions, the couplings of the Higgs bosons 
which involve heavy and light chargino/neutralino states are maximal, while the couplings 
involving only heavy or light ino states are suppressed by powers of M 2 / ' \i for \\i\ ^> M 2 or 
powers of /i/M 2 for \\i\ <C M 2 . To the contrary, the gauge boson couplings to charginos and 
neutralinos are important only for higgsino-like states. Thus, in principle, the (higgsino 
or gaugino-like) heavier chargino and neutralinos X2 an d X34) will dominantly decay, if 
phase space allowed, into Higgs bosons and the lighter x states. 

However, in the asymptotic limit where the heavier chargino and neutralino masses 
are very large, m Xi 3> rn Xj , M Hk , M v , the decay widths into Higgs bosons grow as m xv 

T{ Xl - X jH k ) ~ \am Xi [(Gf jk ) 2 + {Gf jk f} (26) 

while the decay widths into gauge bosons grow as w? 






r(x, - XjV) ~ ^-f m v f + [Gf lV f] (27) 



This is due to the longitudinal component of the gauge boson propagators which introduce 
extra powers of the x% four-momentum in the decay amplitudes. The suppression of the 
{G^y) 2 squared coupling by powers of (/x/M 2 ) 2 or (M 2 //x) 2 depending on whether we 
are in the gaugino or higgsino region will be compensated by the power m x /M§ from the 
matrix element squared. Therefore, the branching ratios for the decays of heavy x particles 
into lighter ones and Higgs or gauge bosons will have the same order of magnitude. Of 
course, as usual, the charged current decay modes will be more important than the neutral 
current decay modes. 

This is exemplified in the figures. In both higgsino and gaugino regions, Fig. 5 and 
Fig. 6 respectively, the decays of charginos xf ari d neutralinos X3 4 into lighter charginos 
and neutralinos and Higgs bosons are not the dominant ones. Still, decays into H ± bosons 
will have substantial branching fractions of the order of 20 to 30%. 



5 We assume in our analysis that the heavier charginos and neutralinos are coming from the decays 
of heavier squarks, including top squarks which is dealt as a special case. Decays of charginos and 
neutralinos into sleptons, which can be lighter than squarks, are relevant only if the former particles 
are gaugino-like since the higgsino-slepton-lepton couplings are rather tiny, unless tan (3 is very large in 
which case the decays into f's could play a role. However, in the large tan/3 case, H ± particles can be 
produced directly with large cross sections. 
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Figure 5: Branching ratios of heavier chargino and neutralinos into the lighter ones and 
gauge/Higgs bosons as functions of their masses for tan/5 = 10. The charged Higgs boson 
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3.2 Cascade decays of squarks and gluinos into H ± bosons 

3.2.1 The case of squarks 

We first focus on the scenario where gluinos are heavier than the scalar partners of the 
light quarks and in this case, all gluinos produced in the processes pp — > gg and pp — *■ qg 
will decay into squarks and light quarks. The squarks will subsequently decay into their 
partner quarks and all types of charginos or neutralinos. The heavier charginos and 
neutralinos will then decay into gauge bosons and Higgs bosons and the lighter chargino 
and neutralino states. Fig. 7 shows the final cross sections times branching ratio to obtain 
one charged Higgs boson in the final state, o x BR(^ H^). For the total production cross 
sections, we use the same set of parton densities and the same scale as in section 2.1. 

The particularly favorable case that we will discuss here is as follows. We choose a 
common squark mass of trig = 800 GeV and fix tan/3 = 10. We also fix \x = 150 GeV 
and vary the wino mass parameter M 2 in such a way that for M 2 > 300 GeV we have 
higgsino-like lighter chargino and neutralinos. For the gluino mass, we will assume in the 
first step that the gaugino masses are unified at the GUT scale so that rrig = 3M 2 , and 
in the second step, we relax this universality assumption and fix the gluino mass to a 
constant value of rrig = 900 GeV, while still varying M 2 . In the latter scenario, the gluino 
is not allowed kinematically to decay into top squarks if m~ tx ~ m g , but in the universal 
case since m g ~ 3M 2 > m^ + m t , gluino decays into stops have to be taken into account. 
In the analysis, we will treat all squarks on equal footing and assume that the branching 
fraction for top squarks decaying into the heavier charginos and neutralinos are the same 
as for the scalar partners of light squarks. [This is a rather conservative approach since 
we have seen in the previous section that the branching ratios for stop decays into bottom 
quarks and the heavier chargino states can be larger.] The charged Higgs bosons masses 
are chosen to be M H ± = 180, 200 and 300 GeV. 

In the universal scenario of Fig. 7 [upper panel], we see that the cross sections times 
branching ratios for iJ ± final states exceed the level of 0.1 pb in most of the displayed 
M 2 range. This means that with a luminosity of J C = 300 fb" 1 which is expected to be 
collected at the LHC, around ~ 30.000 charged Higgs bosons can be produced through 
the cascade decays in this scenario. For small M 2 values, the states X34 an d xt are n °t 
heavy enough for the decays into H ± bosons to occur, in particular for large M#± . When 
these decays are allowed, a x BR(— * iJ ± ) values of the order of 1 pb for Mu± ~ 180 GeV 
and 0.3 pb for M#± ~ 300 GeV can be reached. For increasing values of M 2 , the gluino 
mass increases and the cross section for associated squark and gluino production, which 
is the largest in this case, as well as for gluino pair production drop and o x BR(— > H^) 
decreases accordingly; at some stage, only the cross section for squark production survives 
[since rrig is fixed]. The decrease of a x BR(^ H^) with increasing M 2 is also due to the 
more suppressed phase space for q — ► q'xfiixl since for large M 2 , m x o,m ± ~ M 2 . For 
even larger M 2 values, M 2 > 650 GeV, the channel xi ~^ H ± X\ opens up, and since the 
phase space is more favorable, because m x o ~ M 2 /2, a x BR(^ H^) increases again. 

In the non-universal case [lower panel of Fig. 7], with a gluino mass fixed to m g = 900 
GeV, the cross sections times branching ratios for H ± final states can be much larger. In 



the low M 2 region, the situation is similar to the previous case in which the gluino was 
rather light and contributed substantially to the production cross sections. However, for 
increasing M 2 values, a gg and a gg stay constant contrary to the universal case, while the 
phase space for the decays of the heavier charginos and neutralinos into charged Higgs 
bosons increases, a x BR(^ H^) can then reach the picobarn level even for relatively 
large charged Higgs boson masses, M H ± ~ 300 GeV, leading to samples containing several 
hundred thousand charged Higgs bosons events with the high luminosity J C — 300 fb -1 . 

3.2.2 The case of gluinos 

The production cross sections times branching ratios for squarks decaying into gluinos, 

with the subsequent three-body decays of the gluinos into the heavier charginos and 

neutralinos which then decay into the lighter \ states and charged Higgs bosons are 

shown in Fig. 8. Again we set tan/3 = 10 and fix the gluino and squark masses to 

rriq = 1.2m g ,m g = 800 GeV; the H ± boson masses are taken to be M^± = 180,200 and 

300 GeV. Then we vary the higgsino parameter [i and take the universal scenario where 

M 2 and M3 are related, as well as the non-universal scenario where M 2 = 150 GeV so 

that the heavier x states are higgsino-like for large enough values of |/x|. In both cases 

the cross sections for squark and gluino production are constant and the variation of 

a x BR(i/ ± ) is only due to the variation of the branching ratios BR(g — *■ X3 4QQiXtQQ') 
and BR( X ° i4 , X ± ^ X ±#T )X o^ ±) . 

In the universal scenario, with m g ~ M 3 ~ 3M 2 , a x BR(iJ ± ) is relatively large for 
small values of fi and M#±, when the gaugino-like heavy x states are light enough for 
the decays xl ~^ xfH T and xf ~^ Xi2-^ ± t° occur. In the mixed region, /i ~ M 2 , the 
mass difference between the heavy and light x states are too small to allow for decays in 
H ± bosons. For large values of /i, a x BR(iJ ± ) increases to reach values of the order of 
~ 0.1 pb for Mh± ~ 200 GeV [in particular when the additional channels X3 ~> Xi H T 
open up] before it drops out because of the gradually closing phase space for the decays 

g^wxl^w'xt- 

In the non-universal case with M 2 = 150 GeV, the region with fi < 200 GeV is cut- 
away since the mass difference between the heavier and lighter x states is too small for 
the decays into H ± bosons to occur. For larger /i, the situation becomes similar to the 
universal case and a x BR(iJ ± ) values of the order of 0.3 pb can be reached for not 
too large H ± masses and favorable regions of parameter space, leading to a sample of a 
hundred thousand charged Higgs boson events for the luminosity expected at the LHC. 

3.2.3 The case of top squarks 

The cross section for top squark production at the LHC can be rather large when m g > 
m g > m i x + m t] i n this case all squarks and gluinos will decay into t x final states. The 
stops will decay into heavier x states, with rates shown in Fig. 3, and the latter will 
possibly decay into lighter x particles and H ± bosons as shown in Figs. 5 and 6. The 
situation will then be almost exactly similar to gluino decays discussed in the previous 
subsection. The reason is that gluino decays into the heavy x particle occurred mainly 
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through the virtual exchange of t\ and the branching ratios were controlled by the iibxf 
and £1^X34 vertices. The only difference with the previous discussion will be due to the 
smaller value of rrif 1 [since now t\ is not virtual anymore] which will suppress the phase 
space for t\ — > bx 2 ,txs4 decays. 

Note that an additional contribution to the cross section will come from the direct 
production of stop quarks gg/qq — > Mi, which for m^ ~ 500 GeV can reach the picobarn 
level. This may substantially increase the number of charged Higgs bosons in the final 
state through the cascade decay t% —> bxt ~~ » ^X? 2 H + for instance. 

4. Direct H + decays of squarks and gluinos 

4.1 Two— body decays of squarks into Higgs bosons 

If the mass splitting between two squarks of the same generation is large enough, as is 
generally the case of the (£, b) iso-doublet, the heavier squark can decay into a lighter 
one plus a gauge boson V = W, Z or a Higgs boson $ = h,H, A, H . The partial decay 
widths are given at the tree-level by [the QCD corrections to these decay modes have also 
been calculated and can be found in Rcf. p5|l: 



r(gi->g}V) = ^^v^^.^) (28) 

m^q'^) = jmfc4**A 1/2 (4>/&) ( 29 ) 



In these equations, the couplings of the Higgs bosons to squarks, g^q'.®, read in the case 
of neutral Higgs bosons: 



9qiq2H 



As w M w 

1 



As w M w 

9q\qiA — ~9qiq\A - 



M 2 z s 2 e q (2I* - Ae q s 2 w ) sin(a + (3) + 2m q c 2 e q (A q r q 2 + 2l\iir\ 
- M 2 z s 28q (2I 3 q - 4e q s 2 w ) cos(« + j3) + 2m q c 2 g q {A q r\ - 2I^r q 2 ) 



/i + A„(tan/3)" 



(30) 



2s w M w 

with the coefficients r\ 2 as [a is a mixing angle in the CP-even Higgs sector of the MSSM, 
and at the tree-level, can be expressed only in terms of Ma and tan /3] 

. sin a t cos a h cos a h sin a ,„ . 

r = — -H ' ^ = — -5 ' r l = R ' r 2 = a ■ ( 31 ) 

sinp snip cosp cosp 

In the case of the charged Higgs boson, the couplings to squarks are given, in terms of 
the squark mixing matrices ffl, by 



2 ,T 



9*&H± 



E M tt <«** (**)« ( 32 ) 



2s w M w ^ V Jik ii" V Jij 
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with the matrix Cgg/H± summarizing the couplings of the H ± bosons to the squark current 
eigenstates and it is given by 

„ r- I ml tan/? + m^/tan/3 — Myy sin 2/3 m& (Abtan(3 + jj) \ , > 

fbH± ~ VI y m t {A t / tan {3 + y) 2m t m b / sin 2/3 ) ( °° J 

For the couplings of squarks to the W and Z gauge bosons, one has 

9**w = -^- ( C9 ^ ~ C6 ^ ) (34) 

In Fig. 9, we display the branching ratio for the decays of a bottom squark into the 
lightest top squark and a charged Higgs boson, b\ — > tiH~, b\ — > i\H + as a function of the 
parameter /i with three values of M 2 = 200, 300 and 400 GeV. We have fixed tan/3 = 10 
and the sbottom, stop and charged Higgs boson masses to m^ = 832 GeV, m^ = 430 
GeV and M H ± = 200(300) GeV in the upper (lower) panel. The other competing decays 
of sbottoms are decays into bottom quarks and neutralinos and top quarks and charginos 
as well as bottom and gluino final states. We have assumed universal gaugino masses so 
that the gluino mass is given by rrig ~ 3M 2 ; the strong interaction decay channel b\ — > bg 
is therefore only open for M 2 = 200 GeV and is dominant in this case. 

As can be seen, for M 2 > 300 GeV [i.e when there is no phase space for the decay 
channel b\ — > bg to occur], BR(&i — > t\H~) can be substantial for large \i values, /j, > 700 
GeV, possibly exceeding the level of 50%. The reason for this feature, besides the fact that 
for fi > 800 GeV, the b\ decays into the heavier chargino and neutralinos are kinematically 
closed, is that the sbottom-stop-if ± coupling is strongly enhanced and becomes larger 
than the sbottom-bottom-gaugino coupling which controls the sbottom decays into the 
lighter chargino and neutralinos. For smaller values of M 2 , as pointed out earlier, the 
decay b\ —>■ bg becomes accessible and would be the dominant decay channel. 

Another possibility would be the decays of the heaviest top squark into a sbottom 
quark and charged Higgs bosons, t 2 — > b 12 H + : since the i 2 — h mass splitting can be 
large, there might be enough phase space for this process to occur. However, the cross 
sections for t 2 from direct production are not large, and if the gluinos are heavier than 
squarks, they will decay into t 2 t final states only less than ~ 10% of the time. This decay 
channel is therefore less favored than the previous one. 

Note that for light top squarks, m^ < m t + m x o and m b + m ±, and light charged 

Higgs bosons M H ± < m^ —m x o, the three body decay t\ — > bH + Xi [which is mediated by 
virtual chargino, top or sbottom exchange] is accessible and would compete with the other 
possible three-body decay mode t\ — ► bW + \\ and the loop induced and flavor changing 
decay t± — > cxI- In some areas of the parameter space, the stop decay into charged Higgs 



bosons can be dominant; see the discussions in Ref. 26 
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Figure 9: The branching ratios for bottom squarks decaying into top squarks and charged 
Higgs bosons as a function of fi for tan/3 = 10 and M 2 = 200,300 and 400 GeV. The 
charged Higgs boson mass is taken to be M#± = 200 and 300 GeV, in the upper and 
lower panels, respectively. The two squark masses are taken to be m~ b = 832 GeV and 



m h = 430 GeV. 
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4.2 Three body decays of gluinos into charged Higgs bosons 

Finally, we discuss the direct decays of gluinos into top squarks, bottom quarks and 
charged Higgs bosons, mediated by virtual top quark or bottom squark exchanges; Fig. 10. 




bi v / 

Y 



\?T 



\ff" 



Figure 10: The Feynman diagrams contributing to the three-body decay g — »■ t\bH . 

The Dalitz density for this decay mode, taking into account all the masses of the final 
state particles, including the bottom quark, is given by: 



dr 



dx\dx 2 



(g - H-bh) 



aa s 
64tt' 



-m- n 



dT t + drr + 2dr,r 



n> 



(35) 



In terms of X\ = 2E H ±/m g , x 2 = 2E b /m g and the reduced masses \ix = rnx/m g , the 
squared t, b contributions and the tb interference are given by 
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Vtx 2 



ytcpb b ig b\ g + y b Sf3a b ig a\ g 
(36) 



The various couplings have been given in the previous sections and the Yukawa couplings 
of top and bottom quarks are given by: y t = m t /(\^2swM w S[3) and y b = m b /(V2swM w cp). 



To obtain the partial decay width, one has to integrate over x\ and x 2 with boundary 
conditions: 



2VA*ff+ < x i < 1 + \Mh+ ~ (v 7 /^ + V^h 



(37) 



< x 2 < s T 



(3* 
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with 



1 (xi - 2)(xx - 1 — fi b + Hh ~ ^h+) ~ >/A 

(39) 



2 1 — X\ + flH+ 

1 (xi -2)(x 1 -l — fi b + m x - y, H +) + V& 



2 1 — x\ + n H -\ 



and 



A = {hh+ ~ X T, 



itolHi -(x 1 -l+(i b + fi tl - Hh+Y 



(40) 



The branching fraction for the three-body decay, BR(g — > tJ)H~ +t\bH + ) is illustrated in 
Fig. 11 as a function of \x for tan/5 = 10. We have chosen squark masses of rriq = m~ b . = 1 
TeV, a gluino mass slightly lower, rrig = 900 GeV, and the lighter stop mass to be 
rajj = 433 GeV; for the charged Higgs boson mass we take three values: M H ± = 190, 230 
and 310 GeV. Thus in this scenario, all squarks [including bottom squarks] will decay 
into almost massless quarks and gluinos and the latter will dominantly decay into the 
lighter top squarks and top quarks. The three-body decays g — > i{bH~ and g — > i\bH + 
have therefore to compete with a strong interaction two-body decay, which has a large 
phase space in this case. This is the reason why the branching ratio hardly exceeds the 
one percent level, which occurs for large \x values when the tbH ± couplings are enhanced. 
Note that the smallness of the branching ratio is also due to the smallness of the tbH + 
coupling for the chosen value of tan/3; for larger or smaller values of tan/5, the branching 
ratio can be significantly larger [e.g. an order of magnitude larger than in Fig. 11 for 
tan/3 = 30]. 

In spite of the small branching ratio, the number of H ± final states can be rather large 
at the LHC in the kinematical configuration shown in Fig. 11. Indeed, for the chosen 
squark and gluino masses, the cross section for gluinos coming from direct production 
and from squark decays is at the level of ~ 5 picobarn. With branching fractions of the 
order of a few percent, this means that the cross section for H production via gluino 
decays can reach the level of a fraction of a picobarn leading, in the favorable case, to a 
sample of a few ten thousand H + events with an integrated luminosity of / £ = 300 fb _1 . 

Finally, let us note that for charged Higgs bosons lighter than top quarks, the decay 
t — > bH + will occur at the two-body level and the branching fraction can be rather large [of 
the order of 10 percent for tan/3 = 10]. In the kinematical configuration discussed often in 
this analysis, i.e. squarks heavier than gluinos which are heavier than m^+rrit, all strongly 
interacting particles will cascade into stop and top quarks, with the latter decaying into 
charged Higgs bosons and bottom quarks. The number of H + boson can be therefore 



substantial; see also the discussion given in |Tj|. These processes will complement [and 



possibly, also compete with] searches of H^ 1 bosons from top quark decays with these 
particles produced directly in pp collisions. 
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Figure 11: The branching ratios for direct decays of gluinos into bottom quarks, top 
squarks and charged Higgs bosons as a function of \x for tan/3 = 10 and M^± = 190, 230 



and 310 GeV. The squark masses are taken to be rrig 
the gluino mass is rrig = 900 GeV. 



1 TeV and m^ 
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5. Conclusions 

We have analyzed the production of the charged Higgs particles of the unconstrained 
Minimal Supersymmetric Standard Model from cascade as well as direct decays of squarks 
and gluinos, which have large production cross sections at the LHC due to their strong 
couplings. 

Squarks and gluinos can decay into the heavier chargino and neutralinos, which sub- 
sequently decay into the lighter chargino and neutralino states and charged Higgs bosons. 
We have shown that the branching fractions for these decay modes can be rather large 
in some regions of the MSSM parameter space, and that the cross sections times the 
branching ratios for H ± production can be significantly large. This would allow the pos- 
sibility to produce large samples of these states at the LHC, in particular in situations 
[iJ ± masses larger than the top quark mass or intermediate values of tan /3, for instance] 
where their production cross sections in the direct processes [associated production with 
bottom quarks or with tb final states] are rather small. 

We have also discussed the possibility of producing charged Higgs bosons either through 
direct decays of the scalar partners of heavy quarks, as is the case in the process 61,2 — > 
H~t\, or in direct three-body decays of gluinos, g — > t\bH~ and t\bH + . In some favorable 
regions of the parameter space, the branching ratios of these decay channels, in particular 
in the former case, can be large enough to allow for the detection of the charged Higgs 
bosons. For gluino decays, the situation becomes even more favorable for light H ± bosons 
when the two-body decays t — >• bH + are kinematically possible, leading to a surplus of 
events compared to the case pp — > it — > H ± + X. 
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Thus, there is an additional source of charged Higgs particles at the LHC with inter- 
esting signals^ in most cases, since the final states involve missing energy, multi-leptons 
[from the decay cascades] and heavy flavors [b and t quarks]. These signals would help 
in detecting these particles in the difficult environment of the LHC. In our parton level 
analysis of these process, we did not attempt to scan the entire MSSM parameter space, 
discuss the final state topologies and the possible backgrounds, etc., although for the 
cross sections and branching ratios, we tried to be rather conservative [for instance by 
not including ^-factors, etc.]. A detailed analysis taking into account all backgrounds, 
selection and detection efficiencies in a realistic way, which is beyond the scope of the 
present paper, is required to assess in which part of the MSSM parameter space these 
final states can be isolated experimentally. 
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